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Abstract 

Film-like and thread-like systems are respectively defined by the prop¬ 
erty that one of the constituting homogenous phases has a constant 
thickness (5) or a constant normal section (of largest chord 5). The 
stick probability function of this phase, in the limit 5 0, naturally 

leads to the definition of the correlation function (CF) of a surface or 
of a curve. This CF fairly approximates the generating stick probabil¬ 
ity function in the range of distances larger than 6. The surface and 
the curve CFs respectively behave as 1/r and as 1/r 2 as r approaches 
to zero. This result implies that the relevant small-angle scattering 
intensities behave as Vs/q 2 or as Vc/q in an intermediate range of 
the scattering vector (g) and as V/q A in the outermost g-range. Sim¬ 
ilarly to V, pre-factors Vs and Vc simply depend on some structural 
parameters. Depending on the scale resolution it may happen that a 
given sample looks thread-like at large scale, film-like at small scale 
and particulate at a finer one. An explicit example is reported. To 
practically illustrate the above results, the surface and the curve CFs 
of some simple geometrical shapes have been explicitly evaluated. In 
particular, the CF of the right circular cylinder is explicitly evaluated. 
Its limits, as the height or the diameter the cylinder approaches zero, 
are shown to coincide with the CFs of a circle and of a linear segment, 
respectively. 

Synopsis: The scattering intensities of thread-like and film-like systems 
respectively behave as Vc/q and Vs / q 2 in an intermediate range of the 
scattering vector q. The Vc and Vs expressions are reported. 
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1 Introduction 


Materials characterized by one film-like or by one thread-like phase are since 
long known. Examples of the first kind are vesicles in solutions or oil-water- 
snrfactant systems. Examples of the second kind are polymers or amyloid 
proto-filaments (Avdeev et al, 2013) in solution. Small-angle scattering 
(SAS) is one the most suitable tools to characterize these materials when 
the film thickness or the thread diameter is of the nanometer order. In fact, 
based on the knowledge that the SAS intensity [/(g)] of a circular cylinder 
decreases, at intermediate scattering vector (g) values, as 1/g 2 or as 1/g de¬ 
pending on whether the cylinder’s radius is much larger or much smaller then 
its height (Porod, 1982; Kirste & Oberthiir, 1982), one looks for the existence 
of a plateau in the plot of g 2 /(g) or of g/(g) versus g to conclude that the 
sample respectively contains (plane) flat particles or (straight) rod-like ones. 
Actually this conclusion has a more general validity since it also applies to 
the cases of curved film-like and thread-like phases. To the authors’ knowl¬ 
edge, a first attempt to get this more general result, starting from the basic 
equations of SAS theory (Guinier & Fournet, 1955; Kostorz, 1979; Feigin & 
Svergun, 1987), was done by Teubner (1990) in the film-like case. 

To make this derivation more detailed as well as to extend it to the thread-like 
case is the aim of this paper according to the following plan. Section 2 re¬ 
ports the basic definitions and results of SAS theory for three phase samples. 
Section 3 considers the case where one of the constituting phases is a film of 
thickness S. One shows how the limit S —> 0 of the stick probability function 
(SPF) of the film-like phase leads to the definition of a surface correlation 
function [ 75 (r)]. This 7 $(r) function fairly approximates the generating SPF 
in the range r > 5 and it behaves as 1/r at small distances. Furthermore, 
the surface CF is explicitly evaluated in the case of a sphere (§ 3.4.1), a circle 
(§3.4.2), a rectangle (§3.4.3) and a cubic surface (appendix B). The circle 
is the limit of a right circular cylinder as the cylinder’s height goes to zero. 
Hence, the corresponding limit of the cylinder CF must coincide with that of 
the circle. To verify this point one needs to know the explicit expression of 
the cylinder CF while the only chord-length distribution is presently known 
(Gille, 2014). This CF calculation is performed in closed form in terms of two 
elliptic integral functions (see appendix A). Section 4 analyzes the case of a 
tread-like phase. In the limit of vanishing thread diameter, the correspond¬ 
ing limit of the relevant SPF leads to the definition of the curve correlation 
function [ 7 c (?")] that behaves as 1/r 2 as r —>■ 0 and fairly approximates the 
considered SPF if r > 5. The curve CFs of a linear segment and a circle are 
explicitly worked out in §4.1.1 and §4.1.2. In §4.1.1 one also shows that the 
linear segment CF coincides with the limit of the cylinder CF as the cylin- 
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der’s radius goes to zero. The behaviors of the scattering intensities relevant 
to film-like and thread-like phases are discussed in § 5. The existence of an 
intermediate g-range where the two intensities respectively behave as Vs/q 2 
and as Vc/q is proved in § 5.1 and § 5.2. Of course, both intensities behave as 
V/q A in the outer g-range. The analytic expressions of pre-factors Vs and Vc 
are also worked out and three simple illustrations are reported. In §5.3 one 
discusses the behavior of the scattering intensity of a right parallelepiped. If 
the sizes of the sides differ more than two order of magnitudes, one finds that 
the intensity at first behaves as 1/g, then as 1/g 2 and finally as 1/g 4 in agree¬ 
ment with the fact that the parallelepiped looks as a thread when observed 
on a large scale, as a film when the observation scale becomes smaller and, 
finally, as a parallelepiped when the observation scale has become fine enough 
to resolve all the involved lengths. Section 6 draws the final conclusions. 


2 Basic definitions and properties of the SPFs 


Consider a statically isotropic sample made up of three homogeneous phases. 
The ith of these occupies the spatial set Vi, of volume V), and has scattering 
density n,. The total sample occupies the union set V = VjUk^UVs of volume 
V . The volume fraction of the ith phase is fa = Vi/V. The correlation 
function of the sample is given by (Ciccariello and Riello, 2007) 


7 ( r ) = ^ 


i =1 


n(i;j, k)fa{ 1 - fa 
2 (V 2 ) 


■r*(r), i^j^k. 


( 1 ) 


Here n{i\j,k ) denotes the scattering contrast between phase % and the re¬ 
maining pair of phases (j, k ). It is defined as 


n(i; j, k) = (m - rij) 2 + (m - n k ) 2 - (rij - n k f . (2) 

(rf) is the mean square scattering density fluctuation. It is equal to Xa<i<j< 3 ( n * — 
rij) 2 fa<f)j or to k)fa( 1 — fa)/2. Finally, according to the definition 


r*(r) 


Pj,i{r) ~ fa 2 

fa(l - fa) 


1,2,3, 


(3) 


Tj(r) denotes the CF of the 2 th phase. It is determined by the only geometry 
of the ith phase because Pi,i(r) is the stick probability function (SPF) relevant 
to the phase pair ( 2 , 2 ). The SPFs were first introduced by Debye et al. (1957) 
[see, also, Peterlin (1965) and Goodisman & Brumberger (1971)], in analogy 
with the Patterson function, according to the definition 

Pij( r ) = 4^7/d^/ 3 P*( ri )Pi( ri +™)d2h, i,j = 1,2,3. (4) 
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Here Pi(r) is the characteristic function of the set V,; (i.e. Pi(r) is equal to 
one or zero depending on whether the tip of r falls inside or outside V*) and 
Cj denotes a unit vector that spans all possible directions. The first integral, 
it being an angular average, ensures the assumed isotropy of the sample. 
The probabilistic meaning of the SPFs defined by © implies (Goodisman & 
Brumberger, 1971) that 

Pi :j (oo) = (picpj and P id ( 0) = (5) 

where is Kronecker’s symbol. Once the above equalities are substituted 
within equations ([I]) and ([3]) one finds that, whatever i, 

7 ( 0 ) = Tj(0) = 1 and 7 ( 00 ) = Tj(oo) = 0. (6) 


Definitions © and © imply that the r-dependence of the sample CF is 
determined by that of the SPFs that will now briefly reviewed. To this aim 
it is observed that the SPFs can also be written as (Ciccariello et al. 1981) 


p i,j( 


1 

4vrV 



dv 2 S 3 (r 1 + ruj 


r 2 ), 


(7) 


where S 3 (-), as specified by the index value, denotes the three-dimensional 
Dirac function. From this expression follows that the behavior of P^r), as 
r —* 0, mainly reflects the geometrical features of surface S t that separates 
the zth phase from the remaining ones. In fact, it results that 


Pi,i(r) ~ ^( r ) = 


r Si r 2 21 ,; 


+ 


+ 


4D 127T V 241/ 


Mi + ©i 


( 8 ) 


where S', denotes the area of surface S t (Porod, 1951), L, is the total length 
of edges C, present on Si and, finally, 21 ,, H, and ©, respectively are the 
angularity (Mering & Tchoubar, 1968; Porod, 1967; Ciccariello, 1984), the 
curvosity (Kirste & Porod, 1962) and the sharpness (Ciccariello & Sobry, 
1995) of the <S, surface. These quantities are defined as follows 


21 ,; = 


>Ci L 


1 + (ir -«,(£)) cot (cq(f)) d£ + ^ 


47T 2 

PP' 


(9) 


and 


Mi — - 
2 


'Si 


3H 2 (r)-KrAr) 


d^ 




1 1,3/ 


( 10 ) 

(ii) 
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with 


li,h = --(^1 + 27 t (cot a,/ sin a*) (cot j h3 + cot 7^) + 

(7r — 7 m ) [cot7 M — 2(cot cq cot a,)/sin '7'*,^)] + (12) 

(7T - 7*7) [cot 7^ - 2 (cot a* cot a e )/ sin 7^)] - 

(tt - 7^) [cot 7^ - 2(cot a* sin 7^)/(sin 7„ sin 7^ sin a,)]). 


In ®, the first contribution (Ciccariello et al, 1981) refers to the edges 
and oti(i) denotes the dihedral angle value at the edge point with curvilin¬ 
ear coordinate l. The second contribution (Ciccariello & Benedetti, 1982) 
arises from possible points of contact between different branches of surface 
S, . These points are indexed by J, and S)j denotes the corresponding value of 
the Hessian of the function resulting from the difference of the two branches. 
Kirste & Porod (1962) obtained equation (TTOjl . Here H( r) and Kq(t) re¬ 
spectively denote the mean and the Gaussian curvatures of surface S t at 
the point (with position vector) r. The two curvatures are related to the 
minimum [i? m (r)] and the maximum [i?jvf( r )] curvature radius of the surface 
by 


H( r) 
Kg{ r) 


1 / 1 1 \ 

2 Vi? m (r) ' Rm(y) ' ’ 
1 

Rm{ r) Rm{ r)’ 


(13) 

(14) 


Finally, equation ([IT]) was obtained by Ciccariello & Sobry (1995). The outer 
summation there present is performed over all the vertex points of Si and 
the inner one over all distinct triple of edges that enter into the /th vertex. 
Further, 7 ^ denotes the angle formed by the it\i and the fth edge while cq is 
the dihedral angle at the *th edge. It is also noted that, whenever a couple 
of facets do not share a common edge but only a vertex, one has to consider 
the virtual edge along which the facets, once they are prolonged, intersect 
each other. This fact explains the presence in of the sign factor (—l) tjJ ^. 

One should refer to page 65 of the above paper for the full definition. This 
paper also shows that the SPF of a phase bounded by a polyhedral surface 
is a 3rd degree polynomial within the innermost range of distances, i. e. it is 
exactly given by equation (jH]) with H = K G = 0. 

A further property of SPFs is the fact that, whenever parts of the interfaces 
are parallel to each other at a relative orthogonal distance S, the second 
derivative of the relevant SPF shows a finite (Wu & Schmidt, 1974) or a 
logarithmic discontinuity (Ciccariello, 1989) as r —y S. This property was 
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thoroughly discussed by Ciccariello (1991) and is practically useful (see, e.g ., 
Melnichenko & Ciccariello, 2012 and references therein). This property ap¬ 
plies to film-like samples and yields an oscillatory contribution decreasing as 
in the range q > 2 tt/8. For simplicity, however, this contribution will not 
explictly be considered in §5.1. 


3 Film-like phase and surface CF 

Assume now that phase 1 have a film-like structure, i.e. Vi is a connected 
or disconnected set delimited by a left (<S] j) and a right (<S 1>r ) surface at a 
relative constant orthogonal distance 8. Moreover, one also assumes that Sij 
and Si tr be smooth (i.e. they have no edges, vertices and contact points) 
and that 8 be small in the sense that it obeys the inequality 

3«R,= {(1/R 2 m +1/R 2 M ))-'/ 2 , (15) 

where R s denotes the mean value of the curvature radii of Si = S\jUS\ tr , the 
surface bounding V\. Denote the surface midway Sij and Si >r by Sf. Then, 
iq, the position vector of a generic point P of Vi, can uniquely be written as 

D = D +£z>(ri), (16) 

where iq denotes the position of the point intersection of Sf with the straight 
line [with direction z>(iq)] that going through P, is normal to Sf. £ is the 
orthogonal distance of P from Sf. In equation ([7j) with i = j = 1, the 
infinitesimal volume element dtq located at rq can be written, up to terms 
o(<5), as [see, e.g., equation (4.7) of Ciccariello (1991)] 

dvi = [1 + fif (r i)\dSi d£, -8/2 < £ < 8/2 , (17) 

where dS\ is the infinitesimal surface element of Sf located at rq. Using 
equations (1151) . HTTT) and dTJ), one concludes that the SPF of film-like phase 1 
takes the form 

i r 5 / 2 r s / 2 r r r 

PiAr) = J d^ 2 J dch dS 1 ^ d^ 2 (l +£iH(f 1 )) x 

(l + H(r 2 ))8 3 (ri + 6^(f i) + rCj - r 2 - &v(r 2 )). (18) 
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3.1 The surface CF definition 


In the range of distances 0 < r < 5 this function is reliably approximated 
by equation (J8J) that in the present case reads 


2 g T rp 

P u (r) ss *p li5 (r) = fa -rrr-4 ' 


3 S 


41/ 48 V 


(( 3H 2 -K G ) Sf ) if 0 < r < 5, (19) 


since the smoothness of the surface and the smallness of 5 makes the surface 
area approximation Si ~ 2Sf accurate. 

In the range r > 5, the smallness of 5 makes it possible to set £i = £ 2 = 0 
within the integrand of (1181) so that the SPF becomes 


S 2 Sf 

Pi, lM « -yP-ls(r) if r>6, 


( 20 ) 


with 


7 s( r ) = 4^7 / / dSl I ( r i + ru — r 2 ) , r > 0. (21) 


/ J JSf JSf 


Function 7 $(r) is named surface correlation function because it is fully deter¬ 
mined by the surface set Sf. Function 7 $(r) was first introduced by Teubner 
(1990). [Note, however, that this author used ^7, instead of 7^, as nor¬ 
malization factor in front of the integral reported in (1211) .] 

Figures [U and [21 show the accuracy achieved in approximating the bulk SPF 
Pi t ir by the 3rd degree r-polynomial (fl9|) in the range 0 < r < S and by 
the surface CF, via equation (j20|) . in the range r > 5 for two simple particle 
shapes: the spherical shell and the cylindrical disk. The accuracy drastically 
improves as S gets smaller. 


3.2 Properties of the surface CF 

From equation (120|) and the fact that Pi.fir) is dimensionless it follows that 
the surface CF has dimension length -1 . The leading behavior of the surface 
CF as r —» 0 was obtained by Teubner (1990), under the assumption that the 
surface is closed [he., with no boundary] and smooth, by the same procedure 
that Kirste & Porod (1962) followed to get their relation. It reads 


7 s(r) ~ 


1 

2 r 



R 


M 


+ o(r), 


( 22 ) 


where the angular brackets denote the average over Sf. This relation shows 
that the surface CF diverges as r —y 0 so that it cannot be normalized to one 







y(r) y(r) 



Figure 1: (a): the black continuous line plots 7v,ss(^), the CF of a spherical 
shell [it is recalled that in the case of a single particle one has </>i = 1 and 
Pii(r) — 7v,ss( r ); the last function being given by fhfUD ] with inner and outer 
radii equal to Ri — R — 1 and R 0 = 1.2, so that 5 = ( R a — Ri ) = 0.2 and 
S/R = 0.25. [For the R definition see equation (1T5T) ] . The continuous red line 
is the plot of S 2 Sf^fs,s P h{ r , R)/y within the range 5 < r < 2 R, r ys,s P h{ r ,R) 
being the CF of a spherical surface of radius R = 1 [see equation (TUP ]. The 
thick long-dash blue curve is the plot of fPgstV). The dotted blue (hardly 
visible) and red curves show the approximation discrepancy since they plots 
c • (jv,ss(r) - Vi,s(r)) and c • (7v, ss (r) - S 2 Sps !Sph {r, R)/V) with c = 10. 
The dash-dotted black, blue and red broken lines respectively are the plots 
10 • rjy )SS (r), 10-r<Pi iS (r) and 10 ■ r S 2 S f j s , sp h(r, R)/V (red) . (b): as in (a) 
with Rj = 1 and 6 = 0.02 and c = 100 (the dash-dotted curves merge in 
the single red one). The two panels makes it evident that the approximation 
gets more and more accurate as S/R decreases. 


at the origin as it happens for standard CFs. 

As r —>■ oo, the leading behavior of the surface CF is 


7 s(r) ~ 


Si 
V ' 


This immediately follows from equations ([5]) and since 


lim t s(r) 

r—>• co 


™ 52 s 


/ 




5 2 S f 


V S 2 6 2 
5 2 S f V 2 


Si 

V 


(23) 
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It is also noted that the surface CF defined by equation (ITT!) generalizes the 
two dimensional (2D) CF defined as (Ciccariello, 2009) 

7 = 2nS / dc7 / dSl j dSl ^ Tl +n7 ~~ r2 )’ ( 24 ) 

where S is a plane set and Co is a unit vector that spans the unit circle. To get 
the relation between the two quantities in the case where the Sf set, present 
in (12Tj) . is plane, one chooses the plane where Sf lies as the z — 0 plane. 
Then, in (l2T]i one writes 

(^(iq + rCo — r 2 )dcu = (^(rq + tCj\\ — r 2 )5i(r cos 9)d(pd(cos9), 

[where a)|[ is the projection of u ) onto the z = 0 plane] and, after integrating 
over 9 , one gets 

7sCO = ~r~Q — / do) [ dS t [ dSi<$ 2 (iq + ru - r 2 ), (25) 

4tt S f r J J Sf J Sf 

where the outermost integral is performed on the unit circle. Comparing fl25|i 
with (124|) and assuming that Sf = S one finds 

2 r 7 s(r) = q 2 (r). (26) 


Ciccariello (2010) showed that, as r —» 0, q 2 (r) « (1 — r L/nS) where L 
denotes the length of the boundary of S. Thus, relation (125]) . allows us to 
generalize equation (1221) to the case where Sf is an open surface, and has 
therefore a boundary of perimeter Lf, in the following way 


7s(r) ~ — 


2 7T S' 


+ ^((y:-fiy) V 0(r) 88 r ^ 0 - (27) 

Finally it is observed that the moments of the surface CF can be expressed 

s f 


f 


in terms of the moments of the surface if 


v 


0. . In fact, putting 


M S) m = / r 2m+2 7 5 (r)dr, 

Jo 


(28) 


from equation (|2Tj) follows 

1 


Mi 


S,m 


,2m. 


AttS 

1 


du / dSi / dS 2 S 3 (ri + r - r 2 ) (29) 

/ J J Sf J s f 


dSi / dS' 2 (r 2 — ri ) 2m . 


^SfJsj JSf 

The last expression is equal to the sum of products of different moments of 
surface Sf [see appendix B of Ciccariello (2014)]. In particular, if m = 0, one 
simply finds 


Ms,o = 

An 


(30) 
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3.3 Integral expressions of the derivatives of the sur¬ 
face CF 


The comparison of (ITT]) with equation (3.1) of Ciccaricllo et al. (1981) shows 
that the surface CF and the second derivative of Pi.i(r) nearly have the same 
mathematical structure since the only difference, aside from the normaliza¬ 
tion factors, lies in the angular contribution (z>(fi) -a;)(z>(f 2 ) -cu) that is only 
present in the P ll {r) expression. Consequently, many of the elaborations 
worked out for this expression apply, mutatis mutandis, to the surface CF. 
In particular, let £(r 1 ,r) denote the spherical surface with center at rj and 
radius r. Then equation (12T1) can be written as 

7s(r) = -7—^— 2 l d5 i l dS l ^(fi + ru - r 2 )dS 2 . (31) 

4tT Sf r 2 J Sf JE(T U r) JSj 


Omitting the overbar, denoting by C^( ri, r) the intersection curve of E(i*i, r) 
with Sf (it is formed by the points of «Si that are at distance r from point iq 
of Sf) and proceeding as in the just mentioned paper, one finds that r ys( r ) 
takes the form 


7s(r) = 


dSi 


4tt Sf r Js f Jc s (n,r) i _ + rCj{i )) • &{£)) 


A£. (32) 


In this expression, £ is the curvilinear coordinate of the curve points (note 
that the curvilinear coordinate definition is such that the distance between 
two points, having curvilinear coordinates t and £ + d£, is equal to d£), 
0 (ti + r<2i (£)) denotes the unit vector orthogonal to (and pointing outwardly 
to) Sf at the point r 2 (f?) = iq + ru(£) and Cj(£) specifies the direction of the 
vector (r 2 (f?) — iq). 

If one puts 

cos@(iq, r,£) = z>(rq + ru(£)) ■ Cj{£), (33) 


and 

p( r i, r) 

equation (|32]) becomes 


/c E (n,r) r2 sin 0 (iq, r, £) 


d£, 


(34) 


7s (r) 


1 

4-nSf 



/c E (n.r) 7,2 sin 0 (iq, r, £) 


&£ 


1 

AirSf 




(35) 


Thus, 7 g(r) is proportional to (p(.,r)) s that denotes the average value of 
p(iq,r) over Sf. The nth derivative of 7 s(r) takes the form 


7s ( " ) W = 


47t Sf 


do) / dSi / dS' 2 (-o) • V 2 ) n h(ri + r£o - r 2 ), (36) 
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which implies that the reduced integral expressions involve an integration, 
along the same Cs(r 1; r) curve, of an integrand that changes with the deriva¬ 
tive order n. To get the integrand expression, proceeding as in Ciccariello 
(1995), one denotes by R(w, v ) the parametric equation of <Si, and by R u (u, v) 
and R v (u, v ) the corresponding derivatives with respect to u and v. One puts 


^ (0) (ri,r 2 (£),r) 


_ rtf) _ 

R (u(£),v(£)) —ri|| 2 sinO(ri, r,£) ’ 


(37) 


[where r 2 (£) = R (u(£),v(£)~) and f(£) denotes the unit vector tangent to 
C s (ri,r) at r 2 (£)] and 


A (n) (r 1 ,r 2 (£),r) 


R . — R ■ 

H'Il S\ v S\ u * ( o\ _ 1 O 

||R U x R.„|| sinO(ri,r ,£) ’ ’ 


(38) 


where the reported partial derivatives are evaluated at the argument values 
shown on the left hand side. In this way, the integral expression of the nth 
derivative of the surface CF becomes 



1 

ATTSf 




d£ 


under the assumption that the involved partial derivatives exist. 


(39) 


3.4 Examples of surface CFs 

The following subsections report the surface CFs relevant to a spherical sur¬ 
face, a circle and a rectangle, while the surface CF of a cubic surface is worked 
out in appendix B. 


3.4.1 The CF of spherical surface 


Consider a spherical shell of inner and outer radius equal to R and R+5 with 
6 < R. The bulk CF was calculated by Glatter (1982) (see also Fedorova & 
Emelyanov, 1977) and reads 


7v,mM 


" 7T (r 3 —6r (2R 2 +2RS+S 2 )+85(3R 2 +3RS+S 2 )) 

IT 

irS 2 (2R+5) 2 
2 rV 

n(12r 2 R 2 -r 4 -16rR 3 +6S 2 (2R+8) 2 ) 

12 rV 

ir(r—2(R+5)) 2 (r+4(R.+8)) 

12 V 

0 


if 0 < r < 5, 

ifS<r<2R, 

if 2 R < r < 2 R + 5, 

if 2R + S < r < 2R + 26, 

if 2R + 25 < r, 

(40) 
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with V = 47r((R + 5) 3 — i? 3 )/3. 

As S —y 0, the spherical shell becomes a spherical surface of radius R. The 
corresponding CF is easily evaluated by (J35]). In fact, Cs(ri,r) is a circle 
of radius R\ — r y/l — r 2 /4R 2 and cos@(ri,r, £) = r/2R. In this way, from 
(j35|) follows that 


7>S,.sp/i(D R) 


i- if 0 < r < 2R, 

2 r 

0 if 2 R < r. 


(41) 


This expression is equal to lim^o V^v^ ss (r) / Sf5 2 as required by equation 
(I2U]I . Figures [T^ anddJo allows one to appreciate the accuracy achievable in 
approximating 7 v iSS (r) by 5 2 Sf^s.sphij, R)/V for the cases (R, 5) = (1, 0 . 2 )u 
and ( 1 , 0 . 02 )u, respectively. 


3.4.2 The surface CF of a circle 


Consider a circular cylinder of radius R and height h. If one lets h go to zero, 
the cylinder shrinks to a circle of radius R. The surface CF of the circle is 
easily evaluated by equation (HMj) . In fact the value of the 2D CF of a circle 
at a given distance r is proportional to the overlapping area of the outset 
circle and the circle shifted by r. In this way, accounting for the appropriate 
proportionality constant, one finds that the surface CF 7 s,crci{ r , R) of a circle 
is 


n fS,crd(j"j R) 


\I 4R 2 — V ^ I 1 T 

p 2 + — arccos ^ 

4 7T 7 r r 2 R 


if 0 < r < 2 R, 
if 2 R < r. 


(42) 


One easily verifies that 


V r 


cyl 


lim 

h —>0 H Sf 


^dsk{f ? R, h) 


ls,crd(r, R ), 


where 7 dsk{ r ,R^h) denotes the cylinder CF that is given by (11081) . Besides, 
the small distance expansion of (1421) yields 


lS,crd{ r i R) ~ 


1 

2 r 


1 

7 r R 


24v tR 3 


0( r 3 ), 


(43) 


in agreement with equation (1271) . 

Similarly to Fig.[E] Fig.[ 2 ] shows the accuracy achieved in approximating 
7 dsk( r , R, h), the cylinder CF, by the relevant 4^1,5( r ) polynomial within the 
range 0 < r < h and by h 2 S'/ 75 iCrc )(r, R)/V cy i if r > h in the case where the 
cylinder has height equal to h = 0.2u and radius equal to R = lu. 
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Figure 2: The shown curves have the same meaning of those reported in 
Fig. |H They do however refer to a circular cylinder of radius R = In 
and height h = 0.2u. The black, blue and red dash-dotted curves are the 
plots of 10 • r "fdskir, R, h) [see equation fllOSj) ]. 10 • r^i i 5 (r) [see Eq. (fT9j) ] 
and 10 • r h 2 SR)/V cy i [see Eq. (1421) ]. The dotted curves plot the 
approximation error multiplied by c = 15. 

3.4.3 The surface CF of a rectangle 

Consider a rectangle of sides a and b with a < b. Its surface CF 7 s,rct{r, a, b) 
is easily obtained using equation (1261) and recalling that the 2D CF of any 
plane polygon has an algebraic form (Ciccariello, 2009). Hence, using the 
result of this paper, one finds that the surface CF of the aforesaid rectangle 

is 


ls,rct(r] a, b) 


nab— 2{a-\-b)r-\-r 2 
2nrab ‘ 

a 2 -\-2br—2b\/r 2 —a 2 ‘—2ab arcsin(a/r) 

2nrab _ ! 

_ a 2 -\-b 2 -\-r 2 —2bV r 2 —a 2 —2a\/r 2 — b 2 _ 

2nrab 

2a6(arccos(6/r)—arcsin(a/r)) 

2nrab ’ 




0 < r < a, 
a < r < b, 

b < r < y/a? + b 2 , 
\Ja 2 + b 2 < r. 


(44) 


4 Thread-like phase and curve CF 

Consider now the case where one of the sample phases (named again phase 
1 ) is formed by threads having the same normal section of maximal chord 25 
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and area a. Assume that 5 be small (in a sense that will be defined later). 
Then, in the limit a —)■ 0, phase 1 shrinks to a curve C that might have 
branching points. The number of these points is assumed to be negligible. 
The curve C can be parameterized as 

R = R (l)=(X(l),Y(l),Z(t)), (45) 

where £ denotes the curvilinear coordinate of the curve point set at R. It 
is also assumed that C is smooth in the sense that R(f?) is continuously dif¬ 
ferentiable up to the third order (included). Then, the curve at each of its 
points is endowed of two curvature radii: R c (£ ) and Rt(£), respectively named 
curvature radius and torsion radius (Smirnov, 1970). To define these, one 
considers the three mutually orthogonal unit vectors defined as follows 

= = and b(£) = t{£) x n(£), (46) 

and respectively named tangent , normal and binormal vectors. Clearly, f {£) 
specifies the direction of the tangent to C at R(£), h(£) is the unit vector 
orthogonal to f (£) such that these two vectors determine the plane osculating 
the curve at R(T). Besides, the orientation of n(£) is such that the curve looks 
concave at the considered point. The expressions of the aforesaid radii are 
as follows 


R c (£) = l/ 


df(£) 

dt 


and R,({) = 1 j(p^p ■ n(fl) • 


(47) 


The curvature radius is the radius of the circumference osculating the curve 
while the torsion radius specifies how the osculating plane rotates as one 
moves along the curve. It follows that the curvature radius is always positive 
while the torsion radius can be positive at some points of the curve and 
negative at others. 

Similarly to equation (TT5|l one defines the mean of the curvature radii of the 
curve as 

R c =((l/R c 2 + 1/Rt 2 )) c ~ 1/2 , (48) 

where the average is performed integrating all over the length of C. (It is 
noted that C generally consists of many disjoint closed and/or open curves.) 
One goes now back to Pi i i(r), the SPF of the thread-like phase, to analyze in 
more detail the limit o —> 0. First, it is observed that a is small if 5 « R c as 
it will hereafter be assumed. Let r be the position vector of a point of phase 
1 and consider the plane that contains this point and is orthogonal to C. 
Denote respectively by R(^) and by a(£) the intersections of the plane with 
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C and with the thread-like phase. One can then write r = R(£) + £j_ where 
£_i_ lies within er(£), the normal section set of the thread at R(T). Function 
(r) can be written as 

Pi,i(r) = -L- [dh [ d%. ± [ di 2 [ d% ± x (49) 

4?r V Jc Ja{h) JC Ja(l 2 ) 

J dch 63 (R(£i) + + ruj — R(-O) — ^2,±) • 

The volume fraction of the thread-like phase is fairly approximated by <p\ = 
aL/V with L equal to the total length of the threads. From equation (|5]) 
follows that 


Pi,i(0) = aL/V, and P M (oo) = a 2 (L/V) 2 . (50) 


Similarly to the case discussed in §3.1, in the range 0 < r < 26, 
fairly approximated by equation (flUj) . This can be written in the 
form 


Pi,i(r) « ^i,c(r) 


cr L rir5L irr 3 L 
2\ r+ 64 VS' 


r < 26, 


P u (r) is 
algebraic 

(51) 


if one approximates the threads by circular cylinders so that S ~ 2tt 6L, 
H fa 4. an d x G ~ 0. 

In the range r > 25, owing to the smallness of 5, one can neglect the de¬ 
pendence on £i ia and £ 2 , 1 . within the Dirac function. Then, the integrals over 
these variables can explicitly be performed and one finds that 


with 


a, i(r) 


a 2 L 

~ 


7 c (r), r > 25, 


(52) 


lei?) = 4xZ / dc7 J Mx J d ^5 3 (R(£i) + ru - R(O)), r > 0. (53) 

One concludes that the r-dependence of the SPF of a thread-like phase is 
described by function 7 c (r) in the outer distance range and by equation fl5Tf 
in the inner one. 

Function 7 c (r) is determined by curve C and will be named curve correlation 
function. The dimensions of 7 c(r) are length - ". 

First one elaborates equation (l53|i so as to convert it into a one dimensional 
integral. To this aim, the equation is written as 


lc(r) = —^ [ dO [ dS [ di 2 <5(R(4) + R - R(4)), (54) 

4vrLr 2 J c J S (R(4),r) Jc 


16 








where R denotes the position vector of the infinitesimal surface element dS of 
S(R(^!), r), the sphere of radius r centered at the infinitesimal curve element 
d£ i set at R(£i). Note that the modulus of R is r. One denotes the function 
defined by the innermost two integrals of (154]) as 


Pc(h,r) = f d S f d4<5(R(-0) + R-R( 4 ))- ( 55 ) 

Je(r(4),0 Jc 

The function can easily be evaluated choosing a Cartesian orthogonal frame 
having the origin O at the position of d£\ so that R(f'i) = 0. The value of the 
integral can be different from zero if and only if E(0, r) intersects C. Assume 
that the intersection points exist and let P denotes one fo these points. Then, 
one chooses the Z axis along OP. The plane containing dS is perpendicular to 
Z while dt 2 forms an angle, denoted by 9, with axis Z , so that dZ = cos(9)d£ 2 - 
Around P, p c {£ i,r) takes the form f dX f dE f (dZ/cos9)5 1 (X)5 1 (Y)5 1 (Z) 
and its value is 1/ cos(0). Recalling that d£ 2 is parallel to t(£ 2 ), with f 
defined by equation P5k ). one gets 


w( 4 ,r) = (R( 4 )-R (£i))/r, 

cos (6(£i ,r)) = 0(£ h r) ■ f(I 2 ), 


and 


pc(l I,r) = 


E 


(56) 

(57) 

(58) 


«i(<i,r) - Tj(f 2 ) ‘-f cos ( 0 ,(e h r)) ’ 

where index i labels all the points of C that are at distance r from the point 
set at R(^i). It is stressed that equation (1551) determines the direction of Cj 
as well as the curvilinear coordinate value £ 2 (denoted by £ 2 ) owing to the 
fact that the vector on the right side of (156|) must be unimodular. By (154)) 
and (1581) one concludes that the curve CF is proportional to the curvilinear 
average of 1/ cos9(£i,r), i.e. 

7 c(r) = -X{ I X . (59) 

47rr z cosp(., r) c 


The leading asymptotic behavior of the above expression as r —» 0 is worked 
out in appendix C and reads 


7c (r) 


27rr 2 


1 <lA> c + o(r 2 ). 


167r N R c 2/c 


0. 


(60) 


This result shows that the curve CF diverges as 1/r 2 as one approaches the 
origin. The leading behavior of the curve CF as r —>■ oo is immediately 
obtained by equations (1521) and (1517b) and reads 


7 c(r) « L/V + 0(l/r), r —* oo. (61) 
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The moments Mc, m of the curve CF, similarly to the case of the surface CF, 
can be expressed in terms of the moments of the curve. In fact 


POO I 

„2ra+2 / n - -L 


M c , m = / r 2m+2 ~f C (r)dr = 


di x I d£ 2 (R(£ 2 ) -R(£i)) 2m (62) 


AnL 


From this relation follows that Mq.q = L/An. 


4.1 Examples of curve CFs 

Two examples of curve CFs are reported in the following subsections. 


4.1.1 The curve CF of a linear segment 

The first refers to a linear segment of length L. Applying definition (G3]), the 
angle defined by (l57j) is equal to zero. Thus, equation (1591) yields 


lC,ls( r i L ) 


i 

2 7T r 2 


0 


9 \ if 0 < r < L, 

2 7r L r ’ 

if L < r. 


(63) 


By letting the radius of a circular cylinder go to zero, the limit of its CF 
[multiplied by V/La 2 as required by Q52D ] must reproduce (ItxSlh One easily 
checks that this property holds true evaluating the aforesaid limit stemming 
from 7 n di{r, R , h), the cylinder CF given by equation (11091) . Besides this CF 
can be approximated by (15TT) if r < 2 R and by {TiR 2 )^c,is{ r i h) if r > 2 R. 
Fig. [3] shows the achieved accuracy in the case R = A u and h = lu. 


4.1.2 The curve CF of a circumference 


The curve CF of a circumference of radius R also can easily be evaluated. One 
puts ^ = Rip. Given r, the associated ip value is p — arccos ({2R 2 — r 2 )/2R 2 ). 
Then one finds that 9(£,r) = 7 t/2 — (n — ip)/2 so that cos 9(£,r) = . 

Finally, from (l59j) . one gets the circumference CF 


7 C,crc(r) 


R 

7r r 2 V 4i? 2 —r' 2 

0 


if 0 < r < 2R, 
if 2 R < r. 


(64) 


This CF behaves as l/2nr 2 as r —» 0 and diverges as l/[2nR 3 ^ 2 (2R — r) 1 / 2 ] 
as r —> 2 R~. 
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y(r) 



Figure 3: The curves have the same meaning as in Fig.O They do however 
refer to a circular cylinder of radius R = T,u and height h = lu. The black, 
blue and red dash-dotted curves are the plots of 300-r 2 7 n ^(r) [see Eq. (11091) ]. 
300 • r 2 fPi jC (r)) [see Eq. ([HI]) ] and 300 • r 2 /i 7 c,i s (r) [see Eq. (1HH|) ]. The dotted 
blue and red curves are the plot of the error, multiplied by 15, in 0 < r < 2 R 
and r > 2 R, respectively. 


5 Scattering intensity behavior 


The scattering intensity I(q) is the Fourier transform (FT) of equation (JT]) 
times V ( rf) and is, therefore, equal to 


4:7tV(t] 2 ) f°° 

I(q ) = - / r sin(gr) 7 (r)dr. 

Q Jo 

This can be recast into the form (Ciccariello& Ricllo, 2007) 

V 3 


(65) 


I(q) = k)(t)i{l - (t)i)t i(q) 


1=1 


with 


~ 4tt 

Ti(q) = — / r sin(gr)Fj(r)dr. 


( 66 ) 


(67) 


9 Jo 

The existence of the integral is ensured by condition (Eb), while equation 
(JHla) implies the following sum-rule 



2 7T 


( 68 ) 
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responsible for the Porod invariant relation 

roo 

Qp= q 2 I(q)(q)dq = 2n 2 (r l 2 ). (69) 

Jo 

Besides, the Ti(q)s obey the condition TRq) > 0 whatever the scattering vec¬ 
tor q because they are the FTs of the r,(r)s that have a convolution structure 
owing to definitions ([3]) and (J4j). However, these conditions are not sufficient, 
on the basis of (]65]h to ensure the required positiveness of I(q) because some 
of the n(i;j, k) s can be negative. Hence, the last quantities have to obey ap¬ 
propriate constraints to ensure the positiveness of J(g). This conclusion, that 
might look at first surprising, is a consequence of the functional density the¬ 
orem that states that the correlation function is uniquely determined by the 
density value and the interaction potential. On the basis of this statement 
it is clear that the assignment of the sample internal geometry, equivalent to 
assigning the r*(g)s, cannot be fully independent of the n, values. 

Each Pj.j(r) is determined by the length distribution of the chords that have 
both ends within the related phase. In most of the cases the distribution is 
not uniform but has a particular shape that naturally defines some length 
values Li,, Lm as the bounds of the ranges relevant to the small and the 
large distances asymptotic behaviors of Pi,i{r) as well as to the positions of 
possible peaks and shoulders. Since T^r) is obtained from Pi^{r) by sutract- 
ing to this quantity (R\ it happens that is certainly positive around the 
origin owing to (j6jr) and that it can be negative for other r values. The 
aforesaid distances L\ < ... < Lm in turns define some particular scattering 
vector values qi {i — 1,...) through the well known relation q\ = 27r/Lj. The 
qiS are the reciprocal space values where one expects a change in the Tj(g) 
behavior. 

One applies now these considerations to the film-like and thread-like systems 
discussed in the previous sections. There it was shown that the Pi.i(r) SPF 
is certainly endowed of two typical lengths: R s [R c \ [see equations (11151) and 
(ITS]) ] and 5 [25] with R s > 5 [R c > 25] and that Pi.i(r) is fairly approximated 
by a third degree polynomial if r < 5[25] and by a surface [a curve] CF if 
r > 5 [r > 25]. ft is convenient to separately discuss the case of film-like 
systems and that of thread-like ones. 

5.1 Intensity behavior in the film-like case 

To begin with one analyzes first the behavior of Ti (q) under the further 
simplifying assumption that one only has two typical lengths,he. L\ = 5 and 
L 2 = R s . Thus, the T 1 (q) behavior is expected to change as one passes from 
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the g-range [0, 2tt/R s \ to [2n/R s , 2n/S} and, finally, to [2n/6, oo]. By ((2]) 
and the results of 5 3 one finds that 


if 0 < r < 5 
7 s(r) if 5 <r, 


where it has been put 


Vi(r) 


(r)-^ 

M 1 - <i>i) ’ 


and 75 (r) 


5^S ns {r)/V-^ 

M 1 - fa) 


The FT of equation tHTill yields 


f i(«) = fpiW + %(i). 


(70) 


(71) 


(72) 


with 


# 1 ( 9 ) = 

47T 

f 

/ r sin(g r)73!(r)dr 

(73) 


q 

Jo 



47 r 

POO 


7 s(q) = 

q . 

/ r sin(gr) 75 (r)dr. 

(74) 


[It is observed that condition (071) ensures that 7s ( r ) * 0 as r —* oo so that 

integral d74|) exists.] The behavior of Ti(g) within the first g-range [0, 27 t/R s ] 
could be obtained by expanding the function around q = 0. One finds 


^1 (q) ~ ^2 (2n + 1)! + ^ 75 , 2 nJ q 2n + o(q 2N ) (75) 

77—0 ' 


with 


and 


Vfy 2n = / r 2 " + 3 iPi(r )dr 


M ls3n = / r 2 n+ 3 7 s,(r)dr. 


(76) 

(77) 


One sees that M 2n [= (Msp i2n + M 7Sj2n )] is the 2 nth moment of Ti(r). Since 
M<$ t 2 n is 0(5 2n+3 ), the value of M 2n is essentially equal to the corresponding 
momentum of 7 <s( r ) and strongly depends on the way 7 <s( r ) approaches zero 
as r — > oo. Since the last behavior generally is not known, on a general 
ground the only existence and positiveness of M 0 is certain. The existence of 
the other higher moments is certain only if one assumes that 75 ( r ) approaches 
zero by an exponential decrease or, more strongly, that the sample is a dilute 
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collection of particles with a finite maximal size and that the inter-particle 
interference may be neglected (Gninier & Fournet, 1955). 

Consider now the second (/-interval [2tt/R s , 2 ti / S\. If q is close to 2i t/R s , 
contribution fp^g) to Ti (q) can still, though less accurately, be approximated 
by the sum (Mp i0 — A%, 2 <Z 2 / 6 +...) [see (1751) ]. To estimate contribution 7 s (q) 
one must determine its asymptotic behavior as q —> 00 . Integrating equation 
(171)1 by parts one finds that the leading term is 


75(9) ~ 


4vr 8 3 Sfy s (8)/V - 8$ 

q 2 0i(i-</>i) 


cos (q8). 


In this expression, the contribution related to (— Scftl) is canceled by the 
corresponding contribution present in ^(g) since ipi(r) also presents the 
constant term Thus one can write 


7 s(q) ~ 


4vr 5 3 S ns {5)/V 
q 2 0 i(l- 0 i) 


cos (q5). 


(78) 


provides one also writes 


$i(?)» 3 (Vi,„-M ft29 2 /6) + ... (79) 

where the primes denote that the moments have been evaluated by C3 
without subtracting (j)\ from fPi(r). In the sub-interval of the considered 
(/-range such that qR s > 2n and q5 < 1, using (1221) and the relation SSf/V ~ 
4>i, from (1751) one gets 

— / \ 27T 8 , , 

K 717' (80) 

Recalling that ML q is 0(5 3 ) one finds that M ^ 0 oc 8 (q8) 2 /q 2 that, compared 
to the above contribution, can be neglected. One concludes that 


fi (q) ~ 


2 n 5 

q 2 1 - </>i 


.. 2vr 1 

lf T, <q< T 


(81) 


Once one considers the third (/-range, i. e. -f- < g, in direct space the sample 
structure is analyzed on a distance scale smaller that 8 . Then, the limit 
5 —y 0 is no longer valid and the large q behavior of T^g) must be determined 
starting from the FT of r^r) defined by (EJ). An integration by parts and 
the use of (J 8 ]) immediately yields the equivalent of the Porod relation 


fi(</) » 


4 7T Sf 

V <pi(l — (pi) </ 4 


if 


27r 

T <? - 


(82) 
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Confining oneself to three phase film-like systems where film-like phase 1 
fully separates phase 2 from phase 3 (so that the last two phases have no 
common interface), phases 2 and 3 are characterized by the only length R s . 
Then, by the same considerations that yielded (l82]h one concludes that the 
asymptotic leading terms of T 2 (q) and T 3 (g) are 


Ufa) 


fo(i») 


i 2ns f 
l V&(!-&) 9 4 ’ 


if q<°t * = 2,3, 
if g < q i = 2,3. 


(83) 


Substituting the above asymptotic behaviors within equation (J65]) it results 
that the asymptotic leading term of the scattering intensity relevant to a 
three-phase film-like system (with no common interface between phases 2 


and 3) is 

f Es ;f 

/ ( ^{| if 

A A 

A 

(84) 

with 

Vs = 7T n(l; 2, 3) 5 2 

Sf, 

(85) 


V =27 T Sf (( 77,1 — 

n 2 ) 2 + (ni - n 3 ) 2 ). 

(86) 


The main conclusion of this analysis is that: if 5 « R s , the log-log plot of 
the scattering intensity of a three-phase film-like system shows a linear behav¬ 
ior with slope -2 at intermediate q-values and a further linear behavior with 
slope ~4 in the outer q-range (of course, provided q m ax , the largest observed 
scattering vector, obeys to q ma x8 » 27t). Hence, the lower bounds q s and qi 
of the two linear behavior ranges yield an estimate of R s and S through the 
relations R s ~ 27t jq s and 6 ~ 2ti jqi. The intersections of the two straight 
lines with the vertical axis set at log(g) = 0 determine the values of log(Ts) 
and log('P). According to (1851) and f[86j) . the resulting ratio Vs/V determines 
6 if the phase contrasts are known and then, either Vs or V can be used to 
determine the him area Sf. 

The assumption that phases 2 and 3 have no common interface is now re¬ 
laxed. Then, equations (1HTP and (j52|) remain unchanged while the second of 
equations (1831) converts into 


f 2 (g) fa 


2 tt(Si j2 + * 5 * 2 , 3 ) 
V <j) 2 ( 1 - </> 2 ) <? 4 ’ 



(87) 


and a similar one for T 3 (g) with 5*1^ + *S*i i3 — Sf. From these relations and 
(j65j) one obtains the general expressions of Vs and V 


= 7Tn(l; 2, 3) S 2 (Si t 2 + *S*i )3 ), (88) 

= 2 7r[S*i i2 (ni - tt, 2 ) 2 + S*i i3 (t7i - 77 3 ) 2 + S*2, 3 (77 2 - 77 3 ) 2 ]. (89) 


Vs 

V 
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Though these two relations, even combined with the Porod invariant one, are 
no longer sufficient to determine the four structural parameters S, S lt2 , Si ,3 
and S 2 , 3 , they yield however useful general bounds on these parameters. 
Two illustrations are now reported. Although they are naive mathemati¬ 
cal models they illustrate the main points of the previous discussion. The 
first illustration deals with the scattering intensity of a spherical shell whose 
( R , ci)parameters take the values: (l,0.02)u and (0.5,0.1)u. The left panel 

SPHERICAL SHELL 

SPHERICAL SHELL 5=0.02, R= 1 & 5=0.1, R= 1/2 



Figure 4: (a) The red and the blue curve are the plots of the FTs of equation 
gQD with (R, S ) respectively equal to (l, 0 . 02 )u and (0.5,0.1)u, while the 
green and the cyan curves are the FTs of (141 H . The blue and the cyan 
intensities have been multiplied by 10 3 for greater clarity, (b) Kratky plots 
of the previous intensities. The dotted curves refer to the FTs of (ITTTl . The 
blue curves have vertically been shifted by 0.3. The agreement improves as 
6 /R gets smaller. [The q units are u" 1 .] 

of Fig. S] shows, in red and blue, the spherical shells’ scattering intensities 
obtained by the FT of (HUjl and, in cyan and green, the corresponding FTs 
of ^ 7 s,s P h{r,R) [see (HU)], equal to 

-^r-1s, sp h(q) = 12vr R 2 5 sin 2 (g R)/[q 2 (3R 2 + 3 R5 + <5 2 )]. (90) 

[This FT has been obtained by integrating over [0, 2 R\, the total support of 
7<s,sp/i( r ’, -R)-] The thin and the thick continuous straight lines are the plots 
of the leading asymptotic terms given by equation (1841) [and (l85jl and (| 86 |) ]. 
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The figure shows that the intensity relevant to the spherical shell [leaving 
aside the oscillations] shows both a q~ 2 and a q~ A behavior and that the in¬ 
tensity of the associated spherical surface well reproduces the first intensity 
throughout the range [0, 2i r/<5] (he. not only within [2ir/R s , 2^/5]). More¬ 
over it results that the smaller the ratio S/2R the wider becomes the range 
where the q ~ 2 behavior occurs. The right panel of Fig. 0] shows the q 2 I(q) 
versus q plots (also known as Kratky plots) of the two considered spherical 
shells (continuous curves) and of their associated spherical surfaces (broken 
curves). The agreement is much better in the case (1,0.02) owing to the 
smaller S/2R value. 

The second illustration deals with the case of two disk-like cylinders char- 


i(q) 


Disk-like Cylinder 
h=0.02, R=1 & h=0.1, R=0.5 


0.01 0.1 


<? I(q) 




Figure 5: (a): The red and the blue curve are the plots of the FTs of equation 
(11081) with ( R , h ) respectively equal to (l,0.02)u and (0.5,0.1)u, while the 
magenta and the purple curves are the associated FTs of (ITT)) . The blue and 
the magenta intensities have been multiplied by 10 3 for greater clarity; (b): 
Kratky plots of the previous intensities (without any scaling factor). The 
broken curves refer to equation (1911) . 


acterized by (h, R) = (0.02, l)u and (/z, R) = (0.1,0.5)u. The corresponding 
intensities are obtained by numerically Fourier transforming equation (1108|) 
and respectively yield the red and the blue continuous curves shown in Fig. 
EK The associated surface-like intensities are the 3D FTs of (1421) multiplied 
by h 2 Sf/V, i.e. 


IfSj 

V 




21X77 


1 — Jo^f^l?) — 


(91) 


Their plots are shown as magenta and purple broken curves. They practically 
superpose to the red and blue curves throughout [0, 2 i7/5\. The black lines 
are the plots of the relevant leading asymptotic behaviors given by (184)) . 
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One again observes that the g-range, where the intensity behavior is 1/g 2 
becomes wider as h/R decreases. In the outer q range, as expected, the 
surface approximation no longer works since the intensity behaves as 1/g 4 
in agreement with Porod’s law. Fig. [5jo allows one to better appreciate how 
the agreement between the particle and the surface-like intensities improves 
throughout 0 < q < 2n/h as h/R decreases. 


5.2 Intensity behavior in the thread-like case 

Though, in most cases, a thread-like phase exists in presence of another single 
phase, hereafter one assumes that it exists in presence of two further phase 
and that the thread-like phase portion that lies on the interface between 
phases 2 and 3 is negligible so as to have L = L 2 + L%, where L 2 and L 3 
respectively denote the total lengths of the threads that fully lie within phases 
2 and 3. It is also assumed that 


5*2,3 >> 5*1,2 and 5*2, 3 >> 5 *i >3 


(92) 


and that the mean of the curvature radii of S* 2)3 is not smaller than R c: defined 
by (USD- The above assumptions amount to say that the threads are in an 
extended configuration and that the mean distance among them is not smaller 
than R c . These properties make it possible to apply the considerations made 
in §5.1 to the thread-like case. Thus, the equation equivalent of (1741) and 
(17Tb) becomes 


7 c(Q) = 


47T 


a 2 L 


r sin(gr) 7 c(r)dr. 


g 0i(l - </> 1 ) V J 25 
By (|SQD, considering the only leading term, one finds 
2a 2 L f sinfaT) , 2 a 2 L 


(93) 


7 c(?) 


VM 1-<!>1 )q J 2qS ^ 


dx 


VfiRl - 4> i) q 12 


n 


- Si(2g<5) , (94) 


where Si(x) is the sine integral function (Gradshteyn & Ryzhik, 1980) that, 
at small x, behaves as —x. One concludes that 


7 c(o) ~ fi(g) 


7T a 2 L 


V 0i(l - 0i) g 


.. 27T 7T 

lf 7T < « < r 


(95) 


In this g-range the Tj(g)s with i — 2, 3 are negligible since they are 0(<5 3 ). In 
the range q > 2i t/5 the behaviors are similar to those found in the film-like 
case. One finds 


fi(g) « 


(5*1,2 + 5*i, 3 ) 

R0i(l-0i)g 4 ’ 


if ^ < g, (96) 
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and 


^3 = 2,3 if 2 S<q. 

V 0i(l - (pi) g 4 -R c 


(97) 


By (T95|) . f|96]) . (JOTJ) and (1651) the asymptotic behavior of the thread-like scat¬ 
tering intensity is 




^ if < q < S, 

n R c 'l rt ’ 


if J < q. 


(98) 


with 


P c = 


7rn(l;2,3)a 2 L 7rn(l; 2, 3)cr</>iy 


(99) 


and V given by (]89]h 

The main consequence of this analysis is that the log-log plot of the scattering 


Needle-like Cylinder 
h=2, R=0.01 & h=l,R=0.05 

Kq) q*i(q) 




Figure 6: The two panels, similarly to Fig. 0 show the log-log plot and the 
ql(q) versus q plots of the intensities relevant to two needle-like cylinders of 
height (h) and radius (R) as specified at the top of panel (a) as well as their 
thread-like approximation given by (11001) . 

intensity of a thread-like sample may show a linear behavior with slope -1 
at intermediate q values and another linear behavior with slope -4 (he. the 
Porod one) at larger qs. The associated constants Vc and V are related to the 
the section and the length of the thread-like phase, the scattering contrasts 
and the interface surface areas as reported in (199|) and ([88]) . The expressions 
immediately convert to the two phase case ones by setting 77.3 = ri 2 and 
£ 2,3 = 0. In this way the Vc expression reduces to that obtained by Kirste 
& Oberthiir (1988) under the more restrictive assumption that the particles 
are rod-like. In particular, in the two phase case, the determination of Vc 
and V allows one to determine the normal section and the total surface areas 
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of the thread-like phase if one knows the contrast because volume fraction 
(f)\ can be obtained by the Porod invariant value [i.e. equation (1591) ]. 

Figure [6] reports an illustration of these results considering the case of a 
needle-like cylinder, i.e. a cylinder with its height h larger than is diameter 
2 R. The scattering intensity is obtained by the numerical integration of 
(11091) . The intensity associated to the limit a —> 0 of the thread like phase 
is obtained, according to ([52]) . multiplying by ha 2 jV the FT of (1551) . In this 
way one finds 


ha 2 

~V~ 


7 c,M 


2a 

hq 2 


cos (hq) + (hq) Si (hq) 



( 100 ) 


where Si (or) is the sine integral function. The expression on the right hand 
side approaches n h R 2 as q —> 0 and behaves as n 2 R 2 /q in the range qh » 1. 
The figure shows that the region where the intensity behaves as 1 /q enlarges 
as the R/h ratio decreases. Actually, as in the previously reported cases, this 
value must be smaller than few percents for the q behavior to be clearly 
observed. 


5.3 The right parallelepiped case 

It may happen that a given sample appears as made up of threads if observed 
on a coarse length scale and of films if observed on a fine one. Then, its 
scattering intensity will show both a q and a q~ 2 behavior, besides the 
q~ A one at very large qs. A dilute monodisperse and statistically isotropic 
sample made up of right parallelepipeds of sizes a x b x c with a « b « 
c, is a paramount example of this phenomenon, as it will now be shown. 
The monodispersity assumption allows one to confine the attention to the 
behaviors of the CF and the form factor of a single parallelepiped. The CF 
of this particle shape was worked out by Gille (1999) and has an analytic 
form. Since one expects that the parallelepiped looks as a linear segment on 
a length scale greater than b (and smaller than c) and as a rectangle on a 
scale greater than a (and smaller than b), from equations (1521) and (12UD it 
follows that 


7 v,pri(r) « 



-lc,is(r, c) 
ls,vct(r] b, c) 


if b < r < c 
if a < r < b, 


( 101 ) 


where 7 v,pri( r ) is the parallelepiped CF reported by Gille (1999), 7 c,is( r ,c) 
the linear segment CF given by (155|) and ^fs,rct( r ] b, c) the rectangle surface 
CF given by equation (144|) . Figure [8 ]t shows the accuracy achieved by (1101 jl 
in approximating the CF of a right parallelepiped with a = l/100u, b = In 











7v,prl & y. 




Figure 7: (a) Behavior of the parallelepiped CF and its thread and him-like 
approximations; (b) behavior of the scattering intensitiy and its leading g” 1 , 
q~ 2 and g -4 ’’asymptotic” terms; (c) behavior of the ql(q) (blue) and g 2 /(g) 
(red) versus q plots of the parallepiped intensity (continuous curves) and of 
the relevant curve and surface intensity approximations (broken curves). 

and c = lOOu by the relevant him and thread CFs. In fact, the continuos red 
curve is the plot of 7 v, P w( r ) and the broken magenta one that of a- n fs,rct(r] b, c) 
throughout the range a < r < b. The continuous blue curve is the plot of 
10 3 -7 v,pri(r) and the broken cyan one that of 10 3 ■ab-'yc,is(^, c) throughout b < 
r < c (this interval has been linearly mapped over [ 0 . 01 , 1 ], and both curves 
are vertically shifted by 0.3). The form factor Iv tP ri(q) °f the considered 
parallelepiped is the red curve shown in hg. 8 b and is numerically obtained 
evaluating the FT of 7 v,pri{f)- The blue, the magenta and the black linear 
segments, which respectively behave as g _1 , g ” 2 and g -4 , have been obtained 
by equations fl98h ). fl84h l and (198b ) [or (jMb)]. They have been drawn within 
the g-ranges: [2n/c,2n/b], [2n/b, 2ir/a\ and [2-jr/a, 6 • lO 4 ?^ 1 ]. One sees that 
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the g -1 and q~ 2 behaviors are fairly obeyed close to the lower bounds of 
the previous two q intervals. In the proximity of the two upper bounds, 
the intensity deviates from the reported two linear behaviors as it appears 
more evident in figure [8^. Here the continuous and broken blue are the 
plots of qIv, P ri(q ) and of abq^c,is(Q,c) [note that, similarly to the previously 
reported cases, this FT has been evaluated integrating over the full support 
of Jc,is( r , c ) ensure its positiveness] within the [ 27 t/c, 27 t/ 6 ] g-range, and 
the continuous and broken red curves are the Kratky plots of Iv,pri{q) and the 
FT of ab'ys,rct(q ; b, c ) (evaluated over [0, \Jb 2 + c 2 ]) within [27 r/b < q < 2n/a] 
(the last interval has been linearly mapped over the former one). 

6 Conclusion 

The reported analysis has shown that the CF of a thread-like or a film-like 
statistically isotropic sample can be approximated by a 3rd degree polynomial 
in an inner distance range and, externally to this, by the associated curve or 
surface CF. The curve and the surface CFs respectively behave as 1/r 2 and 
as 1/r close to the origin. Consequently, the relevant scattering intensities 
respectively behave as Vc/q or as Vs/q 2 in a range of intermediate qs and 
as "P/g 4 in the outer q range. The V and Vs expressions were determined 
by Porod (1951) and by Teubner (1990). The expression of Vc is new. Both 
Vs and Vc coincide with the expressions obtained by Porod (1982) and by 
Kirste & Oberthiir (1982) in the cases of plane lamellae and circular rods. 
On a practical ground, coefficients V and Vc or Vs can easily be determined 
from the linear portions of the intensity log-log plot. Their knowledge de¬ 
termines, in a (nearly) model independent way, the interface area as well as 
the normal thread section area or the him thickness. Furthermore, a numer¬ 
ical check on the physical consistency of the assumed hlm-like or thread-like 
structure is possible because the lower bound of the g-range where the q~ 2 
or g" 1 behavior occurs must be close to the him thickness value or to the 
thread maximal chord value (estimated from the normal section area), re¬ 
spectively determined by Vs and Vc- The model illustrations, reported in 
§3.4.1, 3.4.2 and 4.1.1, suggest that the linear behaviors in the log-log plots 
are observable if the two typical lengths, 6 and R s for the hlm-like case and 
26 and R c for the thread-like one, differ at least by an order of magnitude. 
This fact is confirmed by Figures 5 and 6. The paper analysis also applies to 
three phase systems. In thie case of hlm-like samples involving three phases, 
the knowledge of coefficients Vs, V and Qp only puts some bounds on the 
involved structural parameters since it is not sufficient to uniquely determine 
them. Whenever the scattering intensity is collected over a wider g-range, 


30 



as it happens using also ultra-small scattering equipments, one might study 
samples that behave in a thread-like way in the innermost g-range and in a 
film-like one in the intermediate g-range (see the model illustration reported 
in §5.3). In this case, a trivial extension of the above analysis makes it pos¬ 
sible to determine both the thread section area and the him thickness of the 
analyzed sample. 


Appendix A: the circular cylinder CF 


The chord length distribution (CLD) of a circular right cylinder of radius R 
and height h is since long known (Gille, 2014). It involves the elliptic integral 
functions E(ip, k) and F(ip, k) as well as the complete elliptic integrals E {k) 
and K (k) [we are adopting here Gradsteyn & Ryzhik’s (1980) definitions]. 
Besides, it is known that the CLD takes two different form depending on 
whether the cylinder has a disk-like form (he. h < 2 R) or a needle-like 
one (he. 2 R < h). In the following two subsection one reports the explicit 
expressions of the CF that so far were never written down. 

To this aim, one first puts 


A! 


VI R2 


g = 


<y?3 


A 2 = Vr 2 - h 2 , A 3 = V4 R 2 + h 2 

1 c — r f R 
24:71 rh.R 2 ’ ^ = 2R’ ^ = r A 3 ’ 

<Pi = arcsin£, ip 2 = arcsin G 
A 2 . rA 3 

2 R r 2 h R 


( 102 ) 

(103) 

(104) 

(105) 


G A (r,R,h)=g 3 (r(47rfl 2 (2 A-r) + (r 2 + 2fl 2 )A!)+ (106) 

8 R 2 (r 2 - R 2 )^:) - 16ft R((r 2 + 4fl 2 )E(0 + (r 2 - 4fl 2 )K(J)) 


and 


G c (r, R,h) = Q 12vr R 2 (h 2 + r 2 - R 2 ) - 

(h 2 - 5r 2 - 26R 2 )A 2 A 3 - 24 R 2 (h 2 + r 2 - R 2 )<p 3 
8 hr({r 2 + 4R 2 )E(<p 4 , i) - (r 2 - 4R 2 )F(<p 4 , i) 


(107) 
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The disk case 


Then, integrating twice the CLD expression (see the deposited part), in the 
disk case [i.e. h < 2 R] one finds that the cylinder CF reads 


7 dsk{r,R,h) 


G A (r , R , h) 

Q \l2nh 2 R 2 + 3r(r 2 + 2R 2 )A 1 + 

f|(3r 4 - h 4 + 6 r 2 {h 2 - R 2 ) + 22 h 2 R 2 - 24 i? 4 ) + 

24lR 2 ((r 2 - R 2 )cp i - (h 2 + r 2 - i? 2 )^)- 

16hR((r 2 + 4R 2 )E(<p 2 , f) + (r 2 - 4i? 2 )F(< / 9 2 ,0) 

G c (r,R,h) 


0 


if 0 < r < h, 


if h < r < 2 R, 

if 2 R<r < V4R 2 + h 2 , 

if V4R 2 + h 2 < r. 

(108) 


The needle case 


Similarly, in the needle case [i.e. 2 R < h], the CF is found to be 


'G A (r,R, h) 

4Q 


7 n< u(r,R,h) = < 


3vr R 2 {2hr — R 2 ) 

"^(r 2 + ■ 

G c (r, R, h) 


2 hr((r 2 + 4i? 2 )E(|) 


— 


0 



if 0 < r < 2R, 


if 2 R <r < h, 
if h < r < \/4R 2 + h 2 , 
if \/4R 2 + h 2 < r. 
(109) 


Appendix B: the cubic surface CF 

In deriving equation (l32j) it was implicitly assumed that the given surface 
intersects its translated image along a curve. In the reality it can happen 
that for some translation directions the intersection be a surface subset. An 
illustration of this phenomenon is shown in Figure [S] that refers to a cubic 
surface of area Sf = 6a 2 . The red polygon, shown in Fig. 8a, is the intersec¬ 
tion set of the outset cubic surface with its image resulting by the translation 
of the cubic surface by rCj. This vector is such that its tip point does not 
lie over one of the cube faces. The case where rCj fully lies over one of the 
cube faces leads to an intersection set that is formed by two surfaces that 
respectively have their borders equal to the red thick continuous and broken 
rectangles shown in Fig.8b. When theonly tip of rCj spans one of the cube 
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faces the intersection set is again a surface whose boundary is given by the 
red rectangle shown in Fig. 8c. To get the CF of the cubic surface it is 
necessary to evaluate integral (12T| imposing the aforesaid constraints on rCj. 
This task was explicitly carried out as detailed in the deposited part of this 
paper. Here one simply reports the surface CF expressions relevant to the 
geometrical configurations illustrated in first three panels of Fig. 8. 

They respectively read 


7 s,cs,a(r,a) = < 


2 a 2 

_ 2 _ 

a 


2 Vr 2 —a 2 


1 _ 
a 

_5_ 

2 r 

-~~2 arccos - 

■k r 

5tt- 12 _ 2 j_r_ 

6 7r r a ' 6a 2 


| 2Vr 2 -2a 2 | 


-arcsm 


2(r 2 +5a 2 ) 
3 7r a 2 r 


arcsm 


r 2 — 2 a 2 
2 r 2 —2 a 2 

7 a 3 — 3 ar 2 

(r 2 —a 2 ) 3 / 2 


+ 


if 0 < r < a, 
if a < r < a/ 2 a, 

( 110 ) 


if v^a < r < a/3 a, 
if a/3 a < r, 


7 s,cs,b(r, a) 


and 


' _i_2_ _|_ 

2 r 7r a 2n a 2 

_ 2+7T _ r | 2 y/r 2 —a 2 \ 

2 nr 2n a 2 n ar 

— arcsin(-) 

nr v r 7 

lo 


if 0 < r < a, 

if a < r < -y/2 a, 
if \/2a < r. 


(Ill) 


7<s, cs,c(h a) 



if 0 < r < a, 
if a < r < \/2 a, 


if \/2 a < r < aV3 a, 


0 if \/3 a < r. 

( 112 ) 

Thus, the CF function of a cubic surface is the sum of equations (jllOj) , 
(lllip and (1 112 H . Its hrst moment is equal to 6a 2 /47r, as required by equation 
(j5UP while its plot is shown in panel (d) of Fig. 8 for the case a — 1. The 
discontinuity, present at r = a, arises from the opposite faces that are parallel 
at distance a. 
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Figure 8: If ru is such that the tip of rCj does not. lie on one of the cube’s 
faces, the intersection of the outset and the translated cubic surfaces is made 
up of the continuous and broken thick segments shown in panel (a). If ruj is 
such that 6 is equal to 0 or -| or 7r, so as to have two pairs of faces sliding 
one over the other, the intersection set is formed by two opposite rectangles 
having as borders the continuous and the broken thick red curves shown in 
the (b) panel that explicitly refers to the case 6 = 0. Finally, if rCj is such 
that its tip point spans one of the cube’s face, the intersection set is formed 
by the rectangle shared by the superposing faces, i.e. the rectangle having 
as border the rectangular curve show in red in panel (c). Panel (d )shows the 
final cubic surface CF given by the sum of the aforesaid three contributions. 
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Appendix C: small distance behavior of the 
curve CF 


In order to prove equation (l60l) one proceeds by expanding the parametric 
equation of the curve around £ = 0 which is the curvilinear coordinate of the 
point taken as the origin of the Cartesian frame. The expansion up to terms 
o(£ 3 ) reads 


R(£) « 


e e e 

6 R, 21 6 R c Rt 2M C 


(113) 


One straightforwardly verifies that vectors f, n and b, obtained by applying 
definitions (T46|) to (111311 are mutually orthogonal unit vectors up to terms 
o(£ 3 ). The condition R(f?) • R(£) = r 2 determines the curvilinear abscissa of 
the points that are at distance r from the origin. As r —» 0, the equation is 
easily solved by iteration, putting l r as first step. The solutions are 


-± 


±r 1 + 


24 R r 


+ 0(r 5 ). 


(114) 


Consider the positive solution. One finds that Cj(£ + ) = R (£+)/r and f(£ + ) = 
rlR j ( p |. In this way, by equation (1FT)) , one obtains 


cos 6 + = (2>(£- 1 _) • t(£ + ) ~ 1 

The result for the negative solution is the 
immediately obtained by (J59]). 


8 R c 2 
same. 


■ + o(r 2 ). (US) 

In this way result (l60|l is 
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